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AIRMOO— A GENERAL  PROGRAM  FOR  THE  RAPID  ASSESSMENT 
OF  AIRBORNE  POLLUTANTS 


1 INTRODUCTION 


Background 

The  National  Environmental  Policy  Act  (NEPA)  established  the  re- 
quirement for  all  Federal  agencies  to  assess  the  environmental  impact  of 
their  programs  and  activities.  Department  of  the  Army  ( DA ) guidance1 2’3 
has  firmly  established  the  imoortance  of  cognizance  and  concern  for  the 
air  quality  on  and  adjacent  to  DA  installations.  Monitoring  and  abate- 
ment methodologies  require  detailed  study  and  the  current  state  of  the 
art  for  .conducting  both  activities,  although  often  expensive,  is  ade- 
quate. The  assessment  methodologies  are,  however,  inadequate.  Assess- 
ment tools  which  can  define  within  resource  limits  the  magnitude  of  an- 
ticipated air  pollution  problems  are  required  by  AR  200-1.  The  user  must 
be  capable  of  attaining  "first  cut"  estimates  of  pollution  levels  in  a 
timely  and  efficient  manner,  so  that  the  Army  can  comply  with  environ- 
mental assessment  requirements.  Data  input  requirements  must  be  kept  at 
a minimum  and  the  system  must  provide  all  data  necessary  to  assure  an  ad- 
equate assessment  at  minimum  cost. 

The  Construction  Engineering  Research  Laboratory  (CERL),  through  the 
development  of  the  Environmental  Technical  Information  System  (ETIS), 
has  developed  the  concept  of  general  models  for  the  assessment  of  envi- 
ronmental impact  as  part  of  the  Regional  Impact  Assessment  System. 

These  general  models  will  be  designed  to  aid  the  environmental  assess- 
ment process  and  will  of  necessity  be  less  sophisticated  than  many 
existing  techniques.  The  greater  sophistication,  cost,  and  user  re- 
quirements inherent  to  many  existing  models  are  only  reasonable  if  the 
resultant  output  is  available  in  a timely  manner  and  is  sufficiently 
more  accurate  or  dependable  to  the  extent  that  it  would  impact  decisions 
being  made.  The  general  models  can  aid  in  the  screening  and  analysis  of 
original  alternatives.  The  more  sophisticated  techniques  are  better 
suited  to  the  analysis  of  specific  cases  and  problem  areas. 


1 National  Environmental  Policy  Act  of  1969,  83  Stat  852,  42  USC  4321  et 
„ seq.  (January  1970). 

*"  Environmental  Protection  and  Enhancement  > Army  Regulation  (AR)  200-1 
^ (Department  of  the  Army,  February  1974). 

Handbook  for  Environmental  Impact  Analysis  > DA  Pamphlet  200-1 
(Department  of  the  Army,  April  1975). 
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The  overall  objective  of  this  study  is  to  develop  easily  acces- 
sible, inexpensive  methodologies  to  aid  the  Army  in  preparing  environ- 
mental impact  assessments  of  air  quality. 

This  report  documents  the  design  and  development  of  a prototype 
computer-aided  method  for  assessing  the  dispersion  of  air  pollutants. 

It  outlines  (1)  the  general  design  of  the  system,  (2)  the  theoretical 
basis  for  estimating  pollutant  concentrations  via  the  system,  and  (3)  a 
preliminary  user  manual  for  computer  implementation  of  the  prototype 
system. 


Approach 

Existing  predictive  techniques4  were  reviewed  and  evaluated  for 
suitability  to  environmental  assessment.  Numerical  algorithms  were  se- 
lected to  meet  user  requirements  and  a plan  for  implementing  them  into 
the  air  quality  impact  model  concept  was  developed.  User  input  require- 
ments and  the  format  of  the  model  output  were  investigated  and  document- 
ed. 


Mode  of  Technology  Transfer 

The  AIRMOD  computer  program  will  be  disseminated  in  accordance  with 
procedures  set  forth  in  AR  18-1,  Management  Information  Systems  Poli- 
cxes.  Objectives , Procedures , and  Responsibilities  (Department  of  the 
Army,  March  1976). 
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David  F.  Menicucci,  Air  Quality  Assessment  Model  (AQUAM)  Data  Reduc- 
tion and  Operations  Guide  (Air  Force  Weapons  Laboratory,  Air  Force 
System  Command,  Kirtland  AFB,  December  1975);  UNIMAP  Models  (Environ- 
mental Protection  Agency,  1972  to  1977). 
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2 THE  FUNCTION  OF  GENERAL  AIR  POLLUTION  ASSESSMENT  MODELS 
AND  THEIR  CONSTRAINTS 


The  Function  of  General  Models 


An  air  pollution  impact  assessment  model  should  be  capable  of  de- 
termining one  or  more  of  the  following:  (1)  the  spatial  distribution  of 
pollutant  concentration,  (2)  * ,e  pollutant  levels  at  a point,  and  (3) 
the  maximum  pollutant  concentration  and  its  point  of  occurrence. 

Spatial  concentration  distribution  may  be  represented  by  isopleths, 
which  are  points  of  equal  concentration.  The  resultant  map  indicates 
the  zones  in  which  specific  levels  of  pollutant  concentration  are  antic- 
ipated. This  sort  of  analysis  requires  extensive  computer  time  and  usu- 
ally requires  specialized  hardware  (e.g.,  a plotter  or  graphics  dis- 
play). When  evaluating  maximum  concentration  levels,  the  user  is 
interested  only  in  the  location  and  magnitude  of  the  maximum  impact. 

The  need  to  determine  the  concentration  at  a point  is  the  result  of  a 
localized  concern  for  the  effects  of  a source  on  a geographic  location 
which,  for  one  reason  or  another,  happens  to  be  a particularly  sensitive 
area  (e.g.,  residential  dwellings  or  habitat  for  a sensitive  species). 
The  assessment  of  impact  is  based  upon  the  concentrations  of  pollutants 
estimated  to  occur  at  a specific  point. 


Constraints  on  Model  Development 

Three  primary  user  requirements  considered  in  the  design  of  the  as- 
sessment model  were:  (1)  responsiveness,  (2)  geographic  coverage,  and 
(3)  cost  considerations. 

The  model  must  respond  to  inquiries  rapidly  and  without  sophis- 
ticated data  acquisition,  input,  or  analysis  requirements.  A model  re- 
quiring time-consuming  acquisition  of  project-specific  data  will  limit 
its  use  by  precluding  the  rapid  analysis  of  many  different  project  al- 
ternatives. 

The  model  should  not  be  regionally  specific,  or,  if  so,  it  should 
incorporate  a data  base  sufficient  for  use  in  the  limited  assessment  of 
alternatives  across  the  entire  United  States.  Regionally  specific 
models,  while  adequate  for  their  particular  area,  cannot  provide  the  DA 
user  with  the  necessary  flexibility. 

The  model  must  be  selected  to  insure  efficient  utilization  of  the 
DA  user's  designated  funds.  A model  which  incurs  expense  beyond  desig- 
nated funding  must  provide  an  Increased  volume  and  quality  of  useful  in- 
formation to  justify  its  expense.  Because  complexity  is  often 


correlated  with  higher  expense  and  with  the  requirement  for  greater  ex- 
pertise, the  typical  DA  user  is  limited  to  the  use  of  general  models, 
and  to  sacrificing  some  accuracy  in  the  interest  of  planning  and  budget 
efficiency.  This  constraint  precludes  the  use  of  many  of  the  current, 
complex  air  pollution  impact  assessment  models. 


Types  of  Air  Pollution  Impact  Assessment  Models 


Air  pollution  impact  assessment  models  in  general  can  be  classified 
by  source  type  and  pollutant  type.  These  models  can  be  employed  for  the 
assessment  of  either  short-term  or  long-term  impacts.  Traditional 
short-term  measures  are  steady-state  10-minute,  1-hour,  or  24-hour  con- 
centrations. Long-term  concentrations  are  monthly,  seasonal,  or  annual 
averages  taking  into  account  historical  meteorological  data  for  the  time 
period  under  study. 


Source  Type 

Three  types  of  general  sources  are  recognized:  (1)  the  point 
source,  (2)  the  line  source,  and  (3)  the  area  source. 

Point  source,  being  the  simplest,  is  the  most  common.  In  this 
case,  the  pollutants  are  emitted  from  an  identifiable  stationary  source 
such  as  an  industrial  facility  gas  stack.  The  emissions  are  discharged 
through  the  stack  and  are  diffused  into  the  atmosphere.  In  the  absence 
of  topographic  effects,  the  resultant  pollutant  pattern  theoretically 
forms  a series  of  concentric  circles  about  the  stack  centerline.  Al- 
though this  very  rarely  occurs  in  nature,  the  idealized  concept  is 
useful  for  visualization  purposes. 

Line  source  can  be  defined  as  an  infinite  series  of  point  sources 
aligned  along  one  directional  axis,  such  as  a highway.  Although  in  re- 
ality line  sources,  like  highways,  are  often  curvilinear,  this  concept 
is  workable  if  the  curved  line  is  viewed  as  a series  of  small,  straight- 
line  segnents. 

An  area  source  is  defined  as  a large  number  of  point  sources  con- 
centrated within  a specific  geographic  area,  such  as  the  city  of 
Chicago. 


Pollutant  Type 

Pollutants  can  be  classified  according  to  either  their  chemical 
properties  or  their  physical  characteristics  and  are  considered  either 
conservative  or  nonconservative  based  upon  whether  they  are  removed  by 
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atmospheric  processes.  Nonconservative  pollutants  include  nitrogen 
oxides,  hydrocarbons,  oxidants,  and  particulates.  The  modeling  of  non- 
conservative pollutant  dispersion  is  a complex  phenomenon  that  is  beyond 
the  scope  of  this  report. 

Pollutants  can  be  treated  as  either  gases  or  particulates  depending 
upon  their  particle  size.  Particles  less  than  20  mic/ons  in  diameter 
behave  like  gases  due  to  their  long  residence  times  in  the  atmosphere, 
while  those  greater  than  20  microns  in  diameter  are  treated  as  particu- 
lates because  of  their  appreciable  settling  velocity.  In  the  10  to  20 
micron  range  it  is  necessary  to  make  a case-by-case  determination,  the 
main  considerations  being  the  relative  effects  of  buoyancy  and  gravity. 
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Special  Cases 

In  addition  to  the  above  model  types,  consideration  must  also  be 
given  to  certain  special  conditions  which  can  modify  concentration  esti- 
mates. These  may  be  of  a physical  or  phenomenological  nature. 

Limited  mixing  conditions:  conditions  occur  when  the  natural  dis- 
persion of  pollutants  is  restricted,  trapping  them  below  a stagnant  in- 
version layer.  This  limited  mixing  condition  could  develop  into  a situ- 
ation known  as  fumigation  in  which  the  plume  actually  touches  the 
ground. 

Topography:  while  gentle  slopes  and  small  undulations  will  have 
minimal  effect  on  the  dispersion  pattern  of  pollutants  into  the  atmo- 
sphere, larger  topographic  variations  can  significantly  affect  ground 
level  pollutant  concentrations. 

Sea  breeze  influence:  the  effects  of  a large  body  of  water  adja- 
cent to  the  study  area  upon  wind  patterns  can  be  significant  as  cool  and 
warm  air  masses  interact. 


3 THE  AIR  POLLUTION  DISPERSION  MODEL  (AIRMOD)  — THEORY 


The  preceding  chapter  outlined  the  function  of  general  air  pollut- 
ant impact  assessment  models  and  their  constraints.  The  following  chap- 
ters will  highlight  the  approach  taken  by  CERL  in  the  development  of  the 
Air  Pollutant  Dispersion  Model  (AIRMOD).  AIRMOD  will  fulfill  the  DA 
users'  need  for  a widely  applicable  assessment  tool.  For  complex  situ- 
ations, more  sophisticated  models  may  be  required  to  estimate  the  air 
quality  impact. 


Basic  Equation 

As  shown  in  Figure  1,  the  x-axis  is  defined  to  extend  horizontally 
away  from  the  source  in  the  direction  of  the  wind;  the  y-axis  is  another 
horizontal  line  running  perpendicular  to  the  x-axis;  and  the  z-axis  ex- 
tends vertically  from  the  base  of  the  emitting  stack. 

Ideal  dispersion  of  a pollutant  in  three  dimensions  then  may  be  re- 
presented by  Gaussian  distributions: 

-V2  f -(Z-HJ2  -(z+He)2  ) 


c(x,y,z)  = 


2to/zu 


[Eq  1] 


where:  c = concentration  at  x,  y,  z (y/nr) 

He  = effective  stack  height  (m) 

Q = emission  rate  (g/sec) 
a = lateral  dispersion  coefficient  (m) 

= vertical  dispersion  coefficient  (m) 
u = wind  speed  at  stack  height  (m/sec). 


For  centerline  concentrations  (y  = 0)  at  ground  level  (z  = 0),  this 
reduces  to: 


Inherent  assumptions  in  this  ideal  case  are 
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1. 

2. 


A steady  state  condition  has  been  reached 
The  pollutant  is  conservative. 


In  order  to  evaluate  the  basic  equation,  it  is  first  necessary  to 
determine  the  values  of  the  dispersion  coefficients  (oy  and  az),  the 
wind  speed  at  stack  height  (u),  and  the  "effective"  stack  height  (h). 


The  dispersion  coefficients  (also  known  as  diffusion  factors)  are 
determined  by  the  equations 


a 


y 


o 

z 


[Eq  3] 
[Eq  4] 


where:  A,B,C,D  = American  Society  of  Meteorological  Engineers 
(ASME)  values  shown  in  Table  1 


x = downwind  distance  from  source  (m) 


o = lateral  dispersion  coefficient  (m) 

* 

oz  = vertical  dispersion  coefficient  (m). 

The  wind  speed  (u)  at  the  stack  height  is  estimated  from  a known 
wind  speed  at  a reference  height  according  to  the  wind  profile  equation: 

u = u0(H/Hq)p  [Eq  5] 

where:  Hq  = reference  height  (m) 

H = stack  height  (m) 

p = exponent  from  Table  2 

uQ  = known  wind  speed  at  HQ  (m/sec). 


Table  1 

ASME  Dispersion  Parameter* 


1 

0.40 

0.91 

0.40 

0.91 

2 

0.36 

0.86 

0.33 

0.86 

3 

0.32 

0.78 

0.22 

0.78 

4 

0.31 

0.71 

0.06 

0.71 

5 

0.40 

0.91 

0.40 

0.91 

6 

0.36 

0.86 

0.33 

0.86 

7 

0.32 

0.78 

0.22 

0.78 

I 

Table  2 

Wind  Profile  Parameters* 


* Reoormended  Guide  for  the  Dispersion  of  Airborne  Effluents , 

second  edition  (American  Society  of  Meteorological  Engineers  [ASME], 
1973). 

♦ Demarrais,  G.  A.,  "Wind  Speed  Profiles  at  Brookhaven  National 
Laboratory,"  Journal  of  Applied  Meteorology  , Vol  16  (1969), 

pp  181-189. 
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The  effective  stack  height  is  the  actual  stack  height  plus  the  ver- 
tical distance  the  pollutants  rise  due  to  buoyancy  and  momentum  before 
beginning  to  disperse. 

Momentum  Rise 


Momentum  rise  is  usually  insignificant  unless  the  exit  velocity  is 
greater  than  10  m/sec.  The  empirical  equation  reported  by  Briggs5 
was 


AH  = 1.89D 


[Eq  6] 


where:  D = stack  diameter  (m) 

R = ratio  of  stack  gas  exit  velocity  to  wind  speed  (m/sec) 
H = plume  rise  (m) 

x = downwind  distance  from  source  (m). 


! 

I 

i 
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Buoyant  rise  is  also  determined  by  empirical  equations  published  by 
Briggs.  First,  a buoyancy  flux  (F)  is  calculated  according  to 


,T-T 
2 s a 


where:  D 

% 

v 


2.45vD 


= stack  diameter  (m) 

= ambient  air  temperature  (°K) 

= stack  gas  exit  temperature  (°K) 
= stack  gas  exit  velocity  (m/sec). 


[Eq  7] 


The  equation  for  calculating  rise  will  use  this  buoyancy  flux  term,  but 
the  technique  will  differ  depending  on  the  stability  class  of  the  atmo- 
sphere (see  Appendix  for  a discussion  of  stability  classes  as  they  are 
used  in  AIRMOD). 


Briggs,  Gary  A.,  Some  Recent  Analysis  of  Plume  Rise  Observations 

(Second  Clean  Air  Congress). 
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For  unstable  or  neutral  atmospheres  (stability  classes  1 through 
3),  a critical  distance  (x)  is  calculated  such  that 


x = 14Fff  if  F < 55  -J—  tEq 

sec 


- 4 

x = 34.49F5  if  F > 55  tEq 

secJ 


The  plume  rise  { H ) is  then  determined  by 

1 2 

AH  = *-6F3x7  [Eq  1QJ 

u 

where:  F = buoyancy  flux  (m  /sec  ) 

u = wind  speed  at  stack  height  (m/sec) 
x = downwind  distance  or  3.5  x,  whichever  is  smaller  (m). 

For  stable  atmospheres  (stability  class  4),  the  temperature  gra- 
dient of  the  atmosphere  becomes  important.  Accordingly,  a stability 
factor  (S)  is  defined  by 

S"T^(3F+r)  [Eq  11] 

a 

where:  T = ambient  air  temperature  (°K) 

d 

41-  = atmospheric  temperature  gradient  (°K/m) 
dz 

p = adiabatic  lapse  rate  (0.0098  °K/m). 

For  distance  (x)  less  than  > rise  1s  calculated  as  in  Eq  10. 
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However,  for  larger  distances,  two  estimates  must  be  made 


AH,  = 2 


'1 


AH0  = 


5F 

3 

r8 


9(^ 

1 

4 


[Eq  12] 
[Eq  13] 


The  smaller  of  the 
Now  the  total 


two  (AH^  and  AH2 ) is  selected  for  the  plume  rise. 

effective  height  (Hg)  may  be  computed  as 

H = H + AH  _ . + AH.  . [Eq  14] 

e momentum  buoyant  L H J 


Inversion  Conditions 


When  inversion  conditions  exist,  the  resulting  concentration  dis- 
tribution is  classified  into  three  zones:  (1)  an  unaffected  zone  near 
the  source,  (2)  an  "affected"  zone,  and  (3)  a vertically  uniform  zone 
far  from  the  source.  The  distance  x^^  is  the  value  of  x when 

= 0.47L  [Eq  lb] 

where:  L = height  of  the  base  of  inversion  layer 
x-^  will  be  used  to  define  the  limits  of  each  of  the  three  zones. 

The  unaffected  zone  extends  from  the  source  to  a distance  x,  from 
the  source.  If  the  inversion  layer  is  higher  than  the  effecti ve^stack 
height,  the  pollutant  distribution  in  this  zone  is  unaffected.  Concen- 
trations are  calculated  using  the  basic  equation  exactly  as  if  no  in- 
version existed.  In  actuality,  though,  the  effective  stack  height 
cannot  be  greater  than  the  height  of  the  inversion  layer  since  the  pol- 
lutants would  be  "trapped"  or  confined  to  the  lower  layer.  So,  if  the 
models  for  plume  rise  indicate  a greater  value,  the  height  of  the  in- 
version layer  (L)  is  used  for  effective  stack  height  (H  ) in  the  basic 
equation.  e 

From  x1  to  2x,  is  an  "affected"  zone.  No  model  exists  for  the  pol- 
lutant distribution  in  this  area.  However,  concentrations  are  estimated 
by  assuming  an  exponential  relationship  with  distance  and  extrapolating 
from  the  pollutant  concentration  values  at  Xj^  and  Zx^. 

In  the  vertically  uniform  zone  (beyond  2xj  ) the  inversion  layer  has 
limited  the  vertical  mixing  to  such  an  extent  that  the  pollutant  distri- 
bution has  been  smoothed  out  vertically  and  is  uniform  from  the  ground 
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to  the  inversion  layer.  The  only  dispersion  that  occurs  is  lateral. 
Accordingly,  the  concentration  (c)  in  this  zone  is  given  by 


c(x  ,0,0)  = 


[Eq  16] 


/ZttLuct 


where:  L = height  of  the  inversion  layer  (m) 

Q = emission  rate  (g/sec) 

c = concentration  (g/m3) 

u = wind  speed  at  stack  height  (m/sec) 

= lateral  dispersion  coefficient  (m). 

Topographic  Effects 

Local  topography  can  profoundly  affect  expected  pollutant  concen- 
tration levels  (via  streamlining,  turbulence,  etc.).  Unfortunately,  the 
treatment  of  this  type  of  phenomenon  is  far  beyond  the  scope  of  a gen- 
eral, widely  applicable  model.  However,  a corrective  term  may  be  incor- 
porated to  adjust  for  the  difference  in  elevation  between  the  point  of 
pollutant  emission  and  the  points  for  which  pollutant  levels  are  to  be 
determined.  Differences  in  elevation  cause  an  apparent  increase  or 
decrease  in  the  effective  stack  height.  Thus,  the  new  effective  stack 
height  ( Hn ) is 

H = H + (E  - E ) [Eq  17] 

n e s r 

where:  H = effective  stack  height  (m) 

= new  effective  stack  height  (m) 

E"  = ground  elevation  at  the  pollutant  source  (m) 

E*  = ground  elevation  at  the  pollutant  receptor  (m). 

Long-Term  Averages 

The  basic  equation,  as  previously  noted,  is  applicable  for  short- 
term estimates  where  the  meteorological  conditions  are  assumed  to  be 
constant.  Over  long  periods  this  is  obviously  not  true.  A long-term 
estimate  may  be  obtained,  however,  if  site-specific  meteorological  data 
giving  the  observed  frequency  distributions  for  stability  class  and  wind 
speed  class  are  available.  The  "wind  rose"  data  of  the  National  Cli- 
matic Center  provide  this  information  for  the  16  cardinal  wind 
directions.  Using  these  data,  a long-term  concentration  may  be  calcu- 
lated by  taking  an  average,  weighted  by  frequency  of  occurrence,  of  all 


stability  classes  and  wind  speed  classes.  Mathematically  this  is 


where:  Q = 

f(9,  s,  w)  = 


s,w 

(a  ) = 

z's 


uw  = 


c(x,0)  = 


s=4  w=4 

l l 

s=l  w=l 


2Qf(0,s,w) 

(2i')2(a2)suwTf 


2<”z>s 


emission  rate  (g/sec) 

probability  in  stability  class  s 

with  wind  class  w in  direction  6. 

effective  stack  height  for  stability 

class  s and  wind  class  w (m) 

vertical  dispersion  coefficient  for  stability 

class  s (m) 

mean  wind  speed  for  wind  class  w (m/sec). 


[Eq  18] 


Line  Source 


The  basic  equation  for  a point  source  is  extended  to  a line  source 
by  considering  the  line  to  be  an  infinite  number  of  points,  calculating 
the  contribution  from  each  point,  and  adding.  This  amounts  to  integrat- 
ing the  basic  equation  along  the  line  direction.  For  an  infinite  line 
the  result  of  this  integration  is 


O i Q1 

c(x,0,0)  = (^-)2  o^u  • sin(<}>) 


f Eq  19] 


where: 


Ql  = source  strength  per  unit  length 


ui  vv.  j u r cnyui  pci 

= wind  speed  (m/sec) 

= angle  of  wind  with  line  source 
= vertical  dispersion  coefficient  (m) 


This  is  the  basic  equation  for  a line  source.  Topographic  effects  and 
long  term  calculations  are  analogous  to  those  for  a point  source. 

The  source  strength  for  a line  source  corresponds  to  the  emission 
rf'f0[.a  P01"*  source  except  that  it  depends,  obviously,  on  the  length 
or  the  line.  In  the  basic  line  source  equation,  the  source  strength 
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refers  to  the  emission  rate  per  unit  length.  When  the  line  is  a high- 
way, the  source  strength  is 


[Eq  20] 


where:  N = number  of  vehicles  per  sec  (sec  ) 

Qv  = pollutant  mass  per  vehicle  per  meter  (g/sec-m). 

The  calculation  of  the  dispersion  coefficient  ox  needs  to  be  mod- 
ified because:  (1)  the  x-distance  is  not  measured  in  the  wind  direction 
but  at  an  angle,  (2)  the  initial  release  of  pollutant  is  actually  from  a 
moving  point  resulting  in  a pollutant  cloud  with  a finite  size.  The 
equation  is  ★ r 

a = A(x  + I)B  [Eq  21] 


where:  A = coefficient  from  Table  3 

B = coefficient  from  Table  3 
I = value  from  Table  3 
x*  = distance  divided  by  cos  0Q  (m), 


Table  3 

ASME  Dispersion  Parameter* 


units  in  meters 


tabil ity 
Class 


N 

.91 

0.40 

0, 

.86 

0.33 

0, 

.78 

0.22 

0, 

.71 

0.06 

0 

.91 

0.40 

0 

.86 

0.33 

0 

• 

00 

: 

0.22 

1 

.91 

4. 

.86 

6. 

.78 

12. 

.71 

93. 

.91 

4. 

.86 

6. 

.78 

12. 

* Recommended  Guide  for  the  Dispersion  of  Airborne  Effluents, 

second  edition  (American  Society  of  Mechanical  Engineers  [ASME],  1973) 


4 THE  AIR  POLLUTION  DISPERSION  MODEL  (AIRMOD)  — USAGE 


To  gain  access  to  AIRMOD,  the  user  must  request  a number,  name,  and 
password  from  CERL  at  217/352-6511,  extension  385  (commercial)  or  958- 
7385  (FTS).  These  codes  will  give  the  user  direct  access  to  the  Envi- 
ronmental Technical  Information  System  (ETIS)  of  which  AIRMOD  is  an  ele- 
ment. 


ETIS  will  begin  automatically  when  the  user  has  successfully  estab- 
lished communication  with  CERL's  computer.  The  first  response  is: 

Welcome  to  CERL's  Environmental  Technical  Information  System 

For  a complete  list  of  ETIS  programs,  type  a carriage  return.  For  imme- 
diate access  to  AIRMOD,  type  5.  Successful  access  to  AIRMOD  is  indi- 
cated by: 


WELCOME  TO  AIRMOD 


Selecting  the  Model  Type 

AIRMOD  will  then  request  a model  type.  Model  options  available  to 
the  user  are  (1)  point,  (2)  line,  (3)  short,  (4)  long,  (5)  limit, 

(6)  detailed,  and  (7)  topographic. 

The  selected  options  should  all  be  typed  on  one  line,  followed  by  a 
carriage  return.  Any  other  words  not  in  the  option  list  may  be  in- 
cluded: “point  short  topographic"  could  alternatively  be  typed  as 
“short-term  point-source  model  with  topographic  data."  To  identify  a 
valid  model  the  user  must  include  one  of  the  options  "point"  or  "line," 
and  one  of  the  options  "short"  or  "long."  "Point"  and  "line"  as  well  as 
"short"  and  "long,"  are  mutually  exclusive;  i.e.,  if  "point"  is  speci- 
fied, do  not  select  "line"  also. 

The  "point"  option  specifies  that  a point-source  model  is  to  be 
used.  Data  which  AIRMOD  will  request  is  shown  in  Table  4. 

The  "line"  option  specifies  that  a line-source  model  is  to  be  used. 
In  the  usual  case,  the  "line"  is  a highway,  so  it  is  assumed  that  there 
is  no  plume  rise  since  pollutants  from  vehicles  are  emitted  at  tem- 
peratures near  the  ambient  temperature  and  with  low  vertical  components 
of  velocity.  The  data  which  AIRMOD  will  request  for  this  model  is  shown 
in  Table  5. 

"Short"  is  the  option  for  short-term  calculations.  No  additional 
data  is  needed  other  than  that  required  for  either  the  "point"  or  "line" 
option.  Short-term  calculations  depend  on  atmospheric  stability  class 


Table  4 

Data  Requirements  for  Point-Source  Calculation 


Data 

Uni  ts 

Emission  rate 

gr/sec 

Stack  height 

m 

Stack  diameter 

m 

Stack  gas  exit  velocity 

m/sec 

Stack  gas  exit  temperature 

degrees  °K 

Ambient  temperature 

degrees  °K 

Maximum  positive  lapse  rate 

oK/m 

Table  5 


Data  Requirements  for 

Line-Source  Calculation 

Data 

Units 

Number  of  vehicles  per  hour 

cars 

Pollutant  vehicle/km 

( gr/car)-km 

Angle  of  wind  with  line  source 

degrees 
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and  wind  speed,  so  the  output  from  this  option  is  a table  showing 
concentration  estimates  for  every  stability  class  and  wind  class. 

"Long"  is  the  option  to  indicate  when  long-term  concentrations  are 
required.  When  the  "long"  option  is  selected,  AIRMOD  types: 

Name  of  meteorological  input  file: 

In  response,  the  user  types  either  the  name  of  an  existing  data 
file,  or  a name  to  be  given  to  the  data  file  about  to  be  created.  If 
the  name  of  an  existing  file  is  typed,  AIRMOD  will  read  the  meteoro- 
logical data  from  that  file.  Otherwise,  AIRMOD  will  create  a file  with 
the  user-assigned  name  and  the  user  will  be  requested  to  type  meteoro- 
logical data.  Each  piece  of  data  is  specifically  requested  by  AIRMOD. 

Long-term  concentrations  depend  on  the  average  wind  speed  and  aver- 
age stability  class.  When  this  option  is  selected,  the  probability  of 
occurrence  of  each  stability  class  and  wind  class  must  be  input.  Input 
can  be  provided  by  (1)  using  the  Unix  "ed"  command6  to  create  a file  of 
data,  or  (2)  inputting  data  interactively  as  AIRMOD  requests  it. 

If  "ed"  is  used  to  create  a file  of  meteorological  data,  it  should 
contain  the  following  information. 

1.  For  each  stability  class  there  should  be  16  lines  of  data. 

Each  line  corresponds  to  a particular  wind  direction.  Line  1 is  N,  line 
2 is  NNE,  line  3 is  NE,  and  so  on  around  the  compass  in  clockwise  order. 

2.  On  each  line  there  must  be  six  numbers.  These  correspond  to 
frequencies  for  a particular  wind  class  1 through  6. 

3.  The  whole  file  then  should  contain  four  groups  of  data.  Each 
group  is  16  lines  long  with  six  numbers  per  line.  These  must  be  kept  in 
order,  i.e.,  stability  class  1 is  first,  direction  0 (north)  comes  first 
within  a stability  class,  and  wind  class  1 is  the  first  entry  on  a line. 

If  "ed"  has  not  been  used  to  prepare  a meteorological  data  file 
ahead  of  time,  AIRMOD  will  prompt  for  each  data  value  as  needed  and  the 
meteorological  file  is  created,  and  will  store  it  under  the  user-as- 
signed name.  These  data  may  then  be  used  during  any  future  AIRMOD  ses- 
sions by  simply  typing  the  user-assigned  name  in  response  to  the  AIRMOD 
query  "Name  of  meteorological  input  file." 


Documents  for  Use  with  the  Unix  Time  Sharing  System  (Western  Electric 
Company,  1975). 


"Limit"  and  "Detailed"  Options 

The  selection  of  the  "limit"  option  will  include  limited  mixing 
conditions  in  the  model  calculations.  Seven  stability  classes  will  be 
used  instead  of  four.  Classes  5,  6,  and  7 will  refer  to  limited  mixing 
classes  (see  Appendix).  This  option  has  no  effect  on  the  long-term  cal- 
culations, since  limited  mixing  conditions  are  not  considered  when  fig- 
uring long-term  averages. 

Normally,  calculations  are  carried  out  for  one  specified  point. 
However,  if  several  points  are  desired,  the  "detailed"  option  is  se- 
lected. To  use  this  option,  an  input  file  containing  the  distances  for 
points  of  calculation  must  be  supplied.  AIRMOO  will  request: 

Name  of  distances  input  file: 

As  for  the  "long"  option,  the  user  then  types  either  (1)  the  name 
of  an  existing  file  containing  distance  data  or  (2)  a name  to  be  given 
to  the  file  about  to  be  created.  AIRMOD  then  either  reads  data  from  the 
named  file  or  interactively  requests  the  necessary  data. 

If  "ed"  is  used  to  create  the  distances  file,  the  format  should  be 
as  follows. 

1.  First  line:  <direction>  <distance>  <optional  topo  data> 

2.  Second  line:  <di recti on>  <distance>  <optional  topo  data> 

3.  Third  line:  <...etc...> 


This  file  should  contain  one  line  for  each  distance  for  which  pol- 
lutant concentrations  are  to  be  calculated.  These  should  be  in  order 
according  to  direction  and  distance  from  smallest  to  largest.  Wind 
directions  are  specified  as  a number  0 to  15  according  to: 


0 - N 4 - E 

1 - NNE  5 - ESE 

2 - NE  6 - SE 

3 - ENE  7 - SSE 


8 - S 12  - W 

9 - SSW  13  - WNW 

10  - SW  14  - NW 

11  -WSW  15  - NNW 


For  example,  a sample  data  file  might  contain: 

0 500 
0 1000 
8 500 
8 1000 


This  file  will  calculate  concentrations  at  distances  of  500  and  1000  m 
to  the  north  (0)  and  to  the  south  (8)  of  the  source.  It  should  be 
noted,  however,  that  short-term  calculations  for  these  two  directions 


****** 


will  De  the  same  unless  topographic  data  are  supplied  Decause  frequency 
data  and  topography  data  provide  the  only  directional  dependence. 


If  "ed"  has  not  been  used  to  prepare  a distance  file  ahead  of  time, 
AIRMOD  will  request  .each  input  individually  and  write  the  file  in  the 
correct  format.  It  may  then  be  referred  to  during  future  AIRMOD  ses- 
sions by  typing  the  user-assigned  name  in  response  to  the  AIRMOD  query 
"Name  of  distances  input  file." 

Topogvccphio 

To  adjust  for  topography,  topographic  data  must  be  input.  If  cal- 
culations are  being  made  for  one  point  only,  AIRMOD  will  request  the 
elevation  at  that  point  and  the  elevation  at  the  source.  If  more  than 
one  point  is  being  used,  however,  topographic  data  for  each  point  is  in- 
cluded with  the  distance  data  discussed  above.  When  AIRMOD  requests: 

Name  of  distances  input  file: 

the  user  then  types  either  (1)  the  name  of  an  existing  file  containing 
distance  data,  or  (2)  the  name  of  the  file  about  to  be  created.  The 
procedure  is  identical  to  that  for  the  "detailed"  option  except  that 
elevation  for  each  point  is  included  in  the  distances  file  after  each 
input  distance.  If  the  "topographic"  option  were  selected,  a sample 
input  file  might  be 


where  the  last  figure  on  each  line  is  the  elevation  (in  meters).  As  for 
the  ’long"  and  "detailed"  option,  AIRMOD  will  request  the  necessary  data 
one  piece  at  a time  if  it  has  not  previously  been  prepared.  Again,  this 
data  may  then  be  referred  to  in  future  AIRMOD  sessions  by  the 
user-assigned  name. 


Thn  HSer  may  ask  for  assi stance  at  any  time  in  an  AIRMOD  session  by 
typing  "?".  For  example,  the  user  may  request  help  on  the  emission  rate 
of  a source  by  typing  "?".  AIRMOD  will  then  suggest  some  reasonable 
emission  rates  for  various  types  of  pollutant  sources.  When  puzzled  or 
in  need  of  further  explanation,  simply  type  "?". 


* t 


X*T 


* 


* 

I 


rtfofc) 


Welcome  to  CERL's: 


Environmental  Technical  Information  System 
What  program?  (type  <cr>  to  see  list)  5 


WELCOME  TO  AIRMOD 


Model  type : short-term  point-source  model 

Distance  from  source  to  point  (meters):  500 

Emission  rate  (gram/ sec):  100 

Stack  height  (meter):  25 

Stack  diameter  (meter):  2 

Stack  gas  exit  velocity  (meter/ sec):  5 

Stack  gas  at  temperature  (degrees  K):  400 

Ambient  temperature  (deg  K):  300 

Max  positive  lapse  rate  (deg/mTT-  ? 

The  maximum  positive  lapse  rate  is  the  temperature  gradient  in  the  atmos- 
phere. This  number  is  not  critical  --  it  is  used  only  for  calculating 
buoyant  plume  rise  in  a stable  atmosphere.  It  is  never  greater  than  0.01 
deg  K.  If  you  type  in  a greater  number,  0.01  will  be  assumed. 

Max  positive  lapse  rate  (deg/m):  J31 

Source  characteristics: 


Emission  rate 

Stack  hei ght 

Stack  diameter 

Stack  gas  exit  velocity 

Stack  gas  temperature 

Air  temperature 

Max  positive  lapse  rate 

DIRECTION  = NORTH  (0) 


100.00  grams/sec 
25.00  meters 

2.00  meters 

5.00  meters/sec 

400.00  degrees  K 

300.00  degrees  K 
0.0100  degrees/meter 


Stab 

Concentration 

in  micro grams/cubic 

cm 

Class 

wind  1 

wind  2 

wind  3 

wind  4 

wind  5 

wind  6 

1 

0.0002 

0.0005 

0.0006 

0.0006 

0.0006 

0.0005 

2 

0.0000 

0.0004 

0.0007 

0.0008 

0.0008 

0.0007 

3 

0.0000 

0.0003 

0.0006 

0.0007 

0.0C07 

0.0006 

4 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

0.0000 

Figure  2.  Transcript  of  AIRMOD  session. 
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5 CONCLUSION 


Documented  in  this  report  are  user  requirements  related  to  air 
quality  impact  analysis  in  quantitative  terms  as  required  for  the  prepa- 
ration of  Environmental  Impact  Assessments  and  Statements.  An  initial 
air  quality  impact  model  design  and  configuration  was  formulated.  This 
formulation  was  based  upon  review  of  existing  models  and  modification  of 
existing  technology  to  meet  Army  users'  needs. 

This  initial  air  quality  impact  model  design  and  configuration  will 
be  studied  further  and  systematized  into  an  operating  model.  Infor- 
mation provided  in  the  report  can  be  used  for  initial  air  quality  impact 
analysis. 
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APPENDIX: 
STABILITY  CLASSES 


When  a gaseous  pollutant  is  released  to  the  air,  it  begins  to  dis- 
perse. The  physical  process  that  transports  the  molecules  of  pollutant 
from  one  point  to  another  is  eddy  diffusion,  or  turbulence.  This  is  a 
random,  or  nearly  random,  motion  of  a large  number  of  molecules  and  is 
physically  related  to  the  temperature  stratification  of  the  atmosphere 
and  to  the  wind  field.  This  leads  to  the  concept  of  stability  classes. 
Stable  layers  are  those  in  which  large-scale  turbulence  is  suppressed, 
while  unstable  layers  allow  relatively  free  mixing. 

Pasquill7  divided  stability  conditions  into  six  classes,  A through 
F (see  Table  Al).  The  determination  of  these  classes  is  shown  in  Table 
Al.  Solar  radiation  (Table  A2)  is  estimated  as  (1)  strong  (if  atmo- 
sphere is  clear  at  solar  altitudes  greater  than  60  degrees),  (2)  mod- 
erate (if  only  solar  altitudes  between  35  and  60  degrees  have  clear 
skies),  and  (3)  slight  (if  only  solar  altitudes  less  than  35  degrees 
have  clear  skies.) 

The  American  Society  of  Meteorological  Engineers  (ASME)8  has  di- 
vided stability  conditions  into  four  classes,  each  one  somewhat  broader 
than  those  defined  by  Pasquill.  These  classes  are  numbered  1 to  4. 

Their  relationship  to  classes  defined  by  Pasquill  is  shown  in  Table  A2. 
For  AIRMOD  purposes,  the  ASME  classes  have  been  used  since,  in  practice, 
the  distinction  between  classes  such  as  "stable"  and  "extremely  stable" 
is  somewhat  arbitrary  and  makes  little  difference  to  the  calculations. 
Three  other  classes  have  been  added  to  allow  AIRMOD  to  incorporate  in- 
version conditions.  These  are  not  true  classes  per  se  but  are  a com- 
bination of  two  of  the  other  classes. 


Inversion  occurs  when  a nonstable  layer  is  overlain  by  a stable 
layer  at  some  height,  known  as  the  mixing  height  (or  lid  height).  Three 
types  of  inversion  are  used  in  AIRMOD:  (1)  very  unstable  overlain  by 
stable,  (2)  unstable  overlain  by  stable,  and  (3)  neutral  overlain  by 
stable.  This  situation  is  depicted  in  Figure  Al. 
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g Pasquill,  F.,  Atmospheric  Diffusion  (Van  Nostrand,  1962). 

Recommended  Guide  for  the  Dispersion  of  Airborne  Effluents  » second 
edition  (ASME,  1973). 
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Table  A1 

Stability  Class  Correlations 


Pasquill  Classes 


ASME  Classes 


Very  unstable 


nstable 


Table  A2 

Estimation  of  Stability  Class 


Surface 

Wind 


2- 3  m/sec 

3- 5  m/sec 
5-6  m/sec 

6 m/sec 


Incoming  Solar  Radiation 


Moderate 


Ni  ght 

4/8  clouds 


3/8  clouds 


Stable  layer 
(Class  4) 

Stable  layer 
(Class  4) 

Very  unstable 
(Class  1) 

Unstable 
(Class  2) 

Class  5 


Class  6 


Stable  layer 
(Class  4) 


Neutral 
(Class  3) 


Class  7 


Hei ght 


Lid  Hei ght 


Figure  Al.  Stability  classes  for  inversion  classes. 

Thus,  AIRMOD  uses  seven  stability  classes  — the  four  defined  by  ASME 
plus  three  pseudo-classes  describing  inversion  conditions.  The  sta- 
bility class  is  tied  closely  to  the  wind  speed.  A given  stability  class 
may  only  exist  within  a certain  range  of  wind  conditions.  So  for  each 
of  the  seven  stability  classes;  six  wind  classes  are  defined  for  each  of 
the  seven  stability  classes  as  shown  in  Table  A3. 


1 Very  unstable 

2 Unstable 

3 Neutral 

4 Stable 

5 Inversion-1 

6 Inversion-2 

7 Inversion-3 
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